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La presente invention concerne une methode permettant de determiner un modele 
de vitesse d' ondes sismiques captees par des recepteurs , couples avec une formation 
souterraine, d'apres des enregistrements multi-deports de ces ondes. 

La methode permet un acces facilite a reformation cinematique avant sommation 
associee aux evenements contenus dans les enregistrements sismiques, qui est combine a 
une technique d r in version cinematique avant sommation en vue de l'imagerie des interfaces 
geologiques du sous-sol. 

Par "information cinematique avant sommation", nous entendons les temps de trajet 
associes aux reflexions enregistrees par des couples source-capteur situes a des distances 
variables Tun de Tautre. A "acces facilite", nous opposons uh pointe manuel des 
evenements sismiques dans les collections multi-deports. Par "methode ^inversion 
cinematique avant sommation", nous entendons une methode permettant, al^partir de 
l'information cinematique extraite des enregistrements sismiques (et non d'une 
approximation de cette information cinematique), de retrouver la geometrie des reflecteurs 
associes ainsi que les vitesses de couche. 

Etat de la technique 

La sismique reflexion est largement utilisee en exploration petroliere, notamment 
pour produire des images de la subsurface a partir de V information contenue dans les ondes 
qui se sont propagees et reflechies sur les discontinuites geologiques du sous-sol. 

Plus precisement, les methodes d'imagerie utilisent 1'information cinematique 
associee aux reflexions sismiques majeures (i.e. les temps de trajet des ondes s'etant 
reflechies sur les discontinuites principales du sous-sol) pour determiner un macro modele 
de vitesse du sous-sol, qui servira lui-meme a transformer les enregistrements sismiques 
temporels en une image profondeur de la subsurface. 
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L'acces a Tinformation cinematique necessaire a la determination du modele de 
vitesse requiert Interpretation des evenements sismiques dans les enregistrements 
sismiques multi-deports. Or les donnees sismiques multi-deports sont generalement 
caracterisees par un mauvais rapport signal sur bruit, d'ou une defaillance des pointes 
5 automatiques classiques et la necessite d'effectuer un pointe des donnees sismiques manuel 
long et couteux. Par ailleurs, pour les campagnes sismiques 3D (predominantes de nos 
jours par rapport aux campagnes 2D), la quantity de donnees a interpreter est un volume 
4D qui peut atteindre des tailles de 1'ordre de la centaine de GigaOctets, voire de 1'ordre du 
TeraOctets, rendant la tache interpretati ve encore plus longue et fastidieuse. 

10 Afin d'eviter cette etape d'interpretation des collections sismiques multi-deports, les 

geophysiciens ont developpe des methodes basees sur des approximations de la geometrie 
des evenements sismiques dans les collections multi-deports. Pour etablir ces 
approximations, ces methodes emettent des hypotheses sur la complexite du sous-sol, 
hypotheses qui peuvent porter sur la geometrie des discontinues geologiques du sous-sol 

15 et sur les variations de vitesse de couches. On citera par exemple la methode decrite dans 
Taner et Kohler (1969) qui suppose un milieu stratifie compose de couches planes 
horizontales a vitesses de couches homogenes et isotropes, celle decrite dans Levin (1971) 
qui etend la methode precedente a des couches pentees. D'autres variantes ont ete 
proposees, mais de fajon generale, les methodes existantes sont basees sur une hypothese 

20 hyperbolique de la courbe des temps de trajet dans les collections multi -deports. Or cette 
hypothese est violee des que les discontinuites geologiques du sous-sol ne sont plus planes 
et/ou des que les vitesses de couche presentent des variations lat6ral.es. 

Nous proposons ici une methode pour approcher au mieux les temps de trajet des 
collections multi-deports dans le cas ou les discontinuites geologiques du sous-sol ont une 

25 geometrie quelconque et ou les vitesses de couche sont moderement variables lateralement, 
et ce, tout en limitant le temps humain requis pour V interpretation des enregistrements 
sismiques. Les temps de trajet obtenus seront alors traites par une methode d'inversion 
cinematique avant sommation qui permet de prendre en compte des cinematiques 
complexes telle que la tomographic de temps de trajet, decrite par exemple dans le 

30 document suivant : 

Bishop, T. et al, 1985, « Tomographic Determination of Velocity and Depth in 
Laterally Varying Media » Geophysics, 50 N°6, 903-923. 
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La methode pourra egalement s'appliquer dans les cas ou les vitesses de couche 
sont variables lateralement, mais de fagon iterative, le modele de vitesse apres mise a jour 
h. une iteration donnee par application d'une methode d'in version cinematique avant 
soramation telle que la tomographic de reflexion precedemment citee servant de modele 
initial pour l'iteration suivante. 

La methode selon l 'invention 

La methode selon l'invention permet de determiner un modele de vitesse d'ondes 
sismiques captees par des recepteurs sismiques couples avec une formation souterraine, en 
reppnse a remission d'ondes sismiques dans le sous-sol par une source sismique, apres 
reflexion sur des interfaces geologiques de la dite formation, a partir d'enregistrements 
multi-deports de ces ondes. 

Pour chaque evenement sismique repere sur les enregistrements et pour chaque 
couche delimitee par les dites interfaces, la methode comporte au moins les etapes 
suivantes : 

a) a partir des enregistrements sismiques avant sommation, on construit une 
collection iso-deport, dont on extrait une information cinematique ou des temps de trajet 
associes a T6venement ; 

*> * 

b) on choisit une gamme de vitesse autour d'une vitesse de reference dans la dite 
couche que Ton echantillonne avec un pas defini ; 

c) pour chaque echantillon de vitesse, on applique une technique d'inversion a 
vitesse fixee afin de determiner, a partir des temps de trajet extraits de la collection iso- 
deport, la geometrie de la dite interface pour l'6chantillon de vitesse en question de fa<jon a 
obtenir une serie de couples interface/ vitesse pour le dit evenement ; 

d) on calcule 1'information cinematique associee a chaque couple interface/vitesse 
obtenu, pour des couples source-recepteur correspondants a des collections multi-deports 
existant dans les enregistrements sismiques ; 

( 

e) pour chaque couple interface/vitesse et pour chaque collection multi-deports 
selectionnee, on evalue la coherence entre les temps de trajet multi-deports ainsi calcules et 
les enregistrements sismiques, et on selectionne pour chaque collection multi-deports la 
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courbe de temps de trajet offrant la coherence maximale avec les enregistrements 
sismiques ; 

f) on applique une methode diversion cinematique avant sommation en utilisant 
les temps de trajet multi-deports obtenus pour l'ensemble des collections multi-deports 
selectionnees, afin de determiner la geometrie et la vitesse de la couche consideree ; et 

g) on itere n fois (n>0) les etapes a) a f) en considerant a chaque iteration le 
modele de vitesse obtenu a l'iteration precedente comme modele de reference pour defmir 
la vitesse de reference de la nouvelle iteration. 

Suivant un mode de mise en oeuvre convenant au cas ou la garnme de vitesse 
choisie ne s'avere pas assez precise a Tissue soit de l'iteration precedente soit des etapes a) 
a f), on realise n fois 1'etape g) avec n>l. 

Suivant un autre mode de mise en oeuvre convenant au cas ou la distribution des 
vitesses varie fortement lateralement et/ou dans la cas ou Ton ne dispose pas d f une 
connaissance a priori suffisamment precise sur la distribution de vitesse dans la couche 
consideree, on realise 1'etape g) sur des gammes de deports de plus en plus larges au fur et 
a mesure des iterations. 

Suivant un autre mode de mise en oeuvre convenant au cas ou la distribution des 
vitesses varie fortement lateralement et/ou dans le cas ou Ton ne dispose pas d'une 
connaissance a priori suffisamment precise sur la distribution de vitesse dans la couche 
consideree, on realise 1'etape g) sur des grilles de collections multi-deports de plus en plus 
fines au fur et a mesure des iterations. 

Suivant un mode de mise en oeuvre convenant au cas ou la geometrie de l'interface 
recherchee est generatrice de triplications, on realise 1'etape e) en considerant les branches 
de triplications dans les collections multi-deports independamment les unes des autres. On 
utilise par exemple des outils de trace de rayons et d'inversion permettant de prendre en 
compte des arrivees multi-valuees. 

Suivant un exemple de mise en oeuvre, a 1'etape a) on construit de preference une 
collection iso-deport a deport nul ou a deport proche. 
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Suivarit un autre exemple de mise en oeuvre, a Tetape c) on applique 
avantageusement une technique d'inversion cinematique a vitesse fix6e telle qu'une 
migration de carte. 

Suivant un autre exemple de mise en ceuvre, on calcule a Tetape d) Tinformation 
cinematique en tra?ant des rayons multi-deports sur Tinterface de chaque couple interface- 
vitesse. 

Suivant un autre exemple de mise en oeuvre, on realise avantageusement Tetape f) 
en appliquant une methode d'inversion cinematique avant sommation telle qu'une 
tomographic de temps de trajet avant sommation. 

Dans sa definition la plus generate correspondant au cas ou le nombre n d'iterations 
a realiser peut eventuelleraent etre nul, la methode s' applique par exemple si la distribution 
de vitesse dans chaque couche est peu variable lateralement et/ou dans le cas ou Ton 
dispose d f un a priori sur la distribution de vitesse dans la couche consid6ree suffisamment 
precis. 

Dans sa definition plus restrictive ou le nombre n d'iterations a realiser est au moins 
egal a 1, la methode s' applique par exemple dans le cas ou la distribution des vitesses dans 
les couches de la formation geologique est fortement variable lateralement et/ou dans le 
cas ou Ton ne dispose pas d'un a priori suffisamment precis sur la distribution cie vitesse 
dans la couche consideree. 

Presentation succincte des figures 

Les caracteristiques et avantages de la methode selon Tinvention, apparaitront plus 
clairement a la lecture de la description ci-apres d'exemples non limitatifs de mise en 
oeuvre, en se referant aux dessins annexes ou : 

- les figures la a lc montrent respectivement une application de la methode a une 
collection multi-deports de type recepteur commun (Fig. la), a une collection recepteur 
commun sur laquelle ont ete supeiposes les temps de trajet multi-deports obtenus avec 
differentes valeurs de vitesse de couche (Fig. lb) et a une collection recepteur commun 
sur laquelle a ete superposee la courbe de temps de trajet multi-deports predite 
montraht le plus de coherence avec Tevenement sismique considere (Fig.lc). 
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Description detaellee 

On dispose d'enxegistrements sismiques que Ton a obtenus au moyen d'un 
dispositif sismique comportant une source sismique emettant des ondes sismiques qui se 
prop agent dans le sous-sol, un ensemble de recepteurs sismiques couples avec le milieu qui 
5 captent les ondes renvoyees par les discontinuity du sous-sol en reponse aux ondes 
emises, et d'un laboratoire d'enregistrement des signaux sismiques captes. 

I - Cas standard 

Nous presentons ci-apres un premier mode de mise en oeuvre de la methode, dans le 
cas ou la distribution des vitesses dans la formation geologique precedant la couche en 
10 question etant connue ou ayant 6t6 estimee par ailleurs, la geometrie de la couche en 
question ne donnant pas lieu a des triplications, la vitesse a determiner etant peu variable 
lateralement ou un a priori suffisamment precis existant sur cette vitesse. On realise les 
etapes suivantes : 

1-1 A partir des enregistrements sismiques avant sommation ou encore multi- 
15 deports, on construit une collection a deport constant, de preference a deport nul (en realite 
une approximation d'une collection a deport nul, connue sous le nom de « stack », 
caracterisee par un meilleur rapport signal sur bruit que les sections iso-deport 
individuelle) ou , a deport proche (correspondant au deport le plus faible des 
enregistrements multi-deports). 

20 Si une interpretation coherente et suffisamment complete de l'evenement sismique 

en cours est realisable dans cette collection, on procede a son pointe directement dans cette 
collection; sinon, afin d'acceder aux temps de trajet associes a Tevenement en question, on 
pourra avoir recours a un detour dans le domaine migre temps ou profondeur suivi d'une 
demigration temps ou profondeur, comme selon la methode decrite dans : 

25 - Ehinger, A., and Lailly, P., 1995, Velocity model determination by the SMART 
method, Part 1: Theory: 65 th Ann. Intemat. Mtg., Soc. Expl. Geophys., Expanded 
Abstracts, pages 739-742. 

1-2 On se donne une vitesse de reference pour la couche en question, ainsi 
qu'une incertitude sur cette vitesse de reference et un pas d'echantillonnage de la 
30 gamme de vitesse ainsi formee; on obtient ainsi une serie d'echanti lions de vitesse pour 
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la couche en question ; la vitesse de reference peut etre quelconque (i.e. variable 
lateralement et verticalement); en particulier, si Ton dispose d'une distribution de 
vitesse obtenue par ailleurs, on pourra utiliser cette vitesse cornme vitesse de reference; 
si Ton ne dispose pas d'information a priori sur la vitesse de la couche, on pourra par 
exemple prendre une vitesse de reference homogene lateralement et verticalement. 

1-3 Pour chacun des 6chantillons de vitesse, on applique une technique, 
d'inversion de type cinematique a vitesse fixee telle qu'une migration de carte par 
exemple, afin de determiner, a partir des temps de trajet extraits de la section a deport 
constant, la geometrie de la dite interface pour l'6chanti lion de vitesse en question ; on 
obtient ainsi une serie de couples interface/ vitesse pour l'evenement en question; on 
trouvera une description de la migration de carte par exemple dans la reference 
suivante : 

- Yilmaz, O., 2001, Seismic Data Analysis - Processing, Inversion and Interpretation 
of Seismic Data: Society of Exploration Geophysicists. 

II peut eventuellement s'averer impossible de trouver un modele d'interface 
expliquant les temps de trajet pour la vitesse fixee en cours (impossibility d'atteindre la 
convergence). Pour une plus grande efficacite de la methode, on peut ne selectionner que 
des couples interface/vitesse pour lesquels l'ecart entre les temps de trajet obsei&es et les 
temps de trajet calcules sur le modele obtenu apres convergence est inferieur a un certain 
seuil. 

A noter que cette inversion h vitesse fixee peut tout aussi bien etre realisee par toute 
autre technique cinematique equivalente. 

1-4 On calcule Tinformation cinematique associee a chaque couple 
interface/vitesse obtenu par exemple par un trace de rayons multi -deports sur Tinterface 
de chaque couple interface/vitesse ainsi obtenu et pour des couples source-recepteur 
correspondant a des collections multi-deports existantes dans les enregistrements 
sismiques ; on choisit de preference de travailler dans des collections point milieu 
commun, en particulier si la vitesse reelle presente des variations laterales non prises en 
compte dans le modele de reference. 
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1-5 Pour chaque couple interface/vitesse et pour chaque collection multi-deports 
selectionnee, on evalue la coherence entre les temps de trajet multi-deports ainsi 
calcules et les enregistrements sismiques ; on peut par exemple employer une technique 
dite de calcul de semblance, pour laquelle on somme 1'energie rencontree dans les 
donnees sismiques le long d r une courbe qui correspondra ici aux temps de trajet multi- 
deports calcules par trace de rayons. Dans un cas ideal, il existera un echantillon de 
vitesse pour lequel la courbe de temps de trajet predite se superposera parfaitement a 
Tevenement sismique considere, ce qui se traduira par une mesure de coherence 
maximale. L est a noter que cette evaluation pourra se faire sur une gamme limitee de 
deport, definie par l'utilisateur ou automatiquement, par exemple par examen de la 
courbe montrant 1'evolution de la coherence en fonction de la gamme de deports. 

Pour chaque collection multi-deports, on selectionne la courbe de temps de 
trajet la plus en accord avec les enregistrements sismiques, et les temps de trajet multi- 
deports constituant cette courbe sont reserves. H est a noter que la courbe de temps de 
trajet offrant la meilleure coherence pourra etre ajustee si necessaire en cherchant dans 
une fenetre verticale predefinie le maximum (ou le minimum en fonction de la polarite 
de I'evenement considere) d'amplitude le plus proche pour chaque trace sismique de la 
collection multi-deport consideree. 

1-6 Les temps de trajet multi-deports ainsi recueillis pour l'ensemble des 
collections multi-deports selectionnees sont alors injectes dans une methode d'inversion 
cinematique avant sommation telle qu'une tomographie de temps de trajet afm de 
determiner la geometrie et la vitesse de la couche en question. 

1-7 II est possible que l'on soit oblige d'iterer les etapes precedentes dans le cas 
ou la coherence avec les enregistrements sismiques n'est pas estimee globalement 
satisfaisante par l'operateur. 

II) Cas particuliers 

* • 

II- 1 Dans le cas ou la distribution des vitesses dans les couches de la formation 
geologique est fortement variable lateralement et/ou dans le cas oil l'on ne dispose pas d'un 
a priori suffisamment precis sur la distribution de vitesse dans la couche consideree, la 
methode selon Tinvention comporte les etapes suivantes : 
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L'ensemble des etapes precedentes est applique de fagon iterative, les temps de 
trajet recoltes au cours d'une iteration etant utilises pour une mise a jour du modele de 
vitesse par une methode d'inversion cinematique avant sommation telle que la 
tomographic de temps de trajet, modele qui sert alors de modele de vitesse d'entree pour 
5 une nouvelle iteration de la methode. 

H est a noter que, afin de s'affranchirau mieux des effets des variations laterales de 
vitesse qui auraient une longueur d'onde plus faible que le deport maximal des collections 
multi-deports, on pourra elargir la gamme de deports consideree par la m£thode au fur et a 
mesure des iterations. Plus precisement, on initiera la methode sur une gamme de deport 

10 limitee, puis la distribution de vitesse trouvee a partir de cette gamme restreinte de deports 
servira de modele de reference pour une nouvelle iteration de la methode lors de laquelle 
on considerera une gamme de deport plus large, et ainsi de suite. H est a noter que le choix 
de la gamme de deports consideree a une iteration donn6e pourra se faire par exemple par 
examen de la courbe montrant revolution de la coherence en fonction de la gamme de 

15 deports. 

Pour les cas ou la distribution des vitesses varie tres fortement lateralement, on 
pourra en plus adopter une approche multi-grilles, plus precisement, on initiera la methode 
sur une grille lache de collections multi-deports (afin de determiner les grandes longueurs 
d'onde des variations de vitesse), puis Ton prendra en compte des grilles de collections 
20 multi-deports de plus en plus fines au fur et a mesure des iterations. 

II-2 Dans le cas ou la geometrie de Tinterface recherchee T est generatrice de 
triplications et si Ton souhaite egalement determiner les zones de la couche consideree 
donnant lieu a ces triplications, il est avant tout necessaire .de disposer de methodes 
d'inversion telles que des methodes d'inversion cinematique (a vitesse fixee et avant 
25 sommation) par exemple permettant de prendre en compte des arrivees multi-valuees. Une 
telle methode est decrite dans la reference suivante : 

- Delprat-Jannaud, F. et Lailly, P., (1995), How to handle multiple arrivals? Journal of 
Geophysical Research, 100, no. Bl, 703-715. 

I 

Les donnees d'entree de ce type de methode sont les temps de trajet multi-values 
30 associes aux evenements sismiques consideres, ainsi que les parametres de rais associes a 
ces temps de trajet. Afin d'obtenir les temps de trajet multi-values et les parametres de rais 
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associes aux donnees a deport constant, on pourra avoir recours a une boucle de migration- 
demigration (comme decrite dans la methode Ehinger and Lailly (1995) precedemment 
citee), en effectuant l'etape de demigration a l'aide d'un trace de rayons permettant le calcul 
d'arrivees multi-valuees. Un tel trace de rayons est decrit par exemple dans : 

- Jurado, R, Lailly, P., and Ehinger, A., (1998), Fast 3D two-point raytracing for 
traveltime tomography: Proceedings of SPIE, Mathematical Methods in Geophysical 
Imaging V, 3453, 70-81. 

Puis, apres avoir constitue les differents couples interface/vitesse par inversion 
cinematique a vitesse fixee sur les donnees temps/parametres de rais ainsi obtenues, on 
calcule Information cinematique, comme decrit dans le cas le plus general, par le biais par 
exemple de traces des rayons multi-deports dans chacun de ces modeles et pour chacune 
des collections multi-deports selectionnees, mais cette fois, adapte aux calculs des arrivees 
multi-valuees tel que celui decrit dans Jurado et al. (1998) precedemment cite. On obtient 
ainsi, pour chaque collection multi-deports s&ectionnee, a la fois des temps de trajet multi- 
values avant sommation et leurs parametres de rais associes. 

Puis, l'analyse de coherence entre les courbes de temps de trajet ainsi predites et les 
evenements consideres dans les donnees sismiques sera realisee tel que decrit dans le cas le 
plus general, mais de facon independante pour chacune des branches des arrivees multi- 
valuees. Les temps de trajet ainsi recoltes pour chacune des branches des arrivees multi- 
valuees seront alors injectes, avec leur parametres de rais associes, dans une methode 
diversion cinematique avant sommation permettant de prendre en compte les arrivees 
multi-valuees, telle que celle decrite dans Delprat-Jannaud et Lailly (1995) precedemment 
citee. 

On a decrit ici une application de la methode a la determination d'un modele de 
Vitesse des ondes sismiques dans une formation souterraine. II est bien evident toutefois 
que la methode peut tout aussi bien s'appliquer a la determination de la vitesse de 
propagation d'autres types d'ondes dans un milieu heterogene presentant des 
discontinuites. ^ 
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REVENDICATIONS 

1) Mdthode pour determiner un modele de vitesse d'ondes sismiques cap tees par 
des recepteurs sismiques couplds avec une formation souterraine, en r6ponse a remission 
d'ondes sismiques dans le sous-sol par une source sismique, apres reflexion sur des 
5 interfaces geologiques de la dite formation, k partir d' enregistrements multi-deports de ces 
ondes, caract&isee en ce qu'elle comporte ay moins les etapes suivantes pour chaque 
evenement sismique repere sur les enregistrements et pour chaque couche delimitee par les 
dites interfaces : 

a) a partir des enregistrements sismiques avant sommation, on construit une 
10 collection iso-deport, dont on extrait une information cinematique ou des temps de trajet 

associes a Tevenement ; 

b) on choisit une gamme de vitesse autour d'une vitesse de reference dans la dite 
couche que Ton echantillonne avec un pas defini ; 

c) pour chaque echantillon de vitesse, on applique une technique d'inversion a 
15 vitesse fixee afin de determiner, a partir des temps de trajet extraits de la collection iso- 
deport, la geometrie de la dite interface pour Techantillon de vitesse en question de fa?on a 
obtenir une serie de couples interface/ vitesse pour le dit evenement ; 

d) on calcule Tinformation cinematique associee a chaque couple interface/vitesse 
obtenu, pour des couples source-recepteur correspondants a des collections multi-deports 

20 existant dans les enregistrements sismiques ; ■ 

e) pour chaque couple interface/vitesse et pour chaque collection multi-deports 
selectionnee, on evalue la coherence entre les temps de trajet multi-deports ainsi calcules et 
les enregistrements sismiques, et on selectionne pour chaque collection multi-deports la 
courbe de temps de trajet offrant la coherence maximale avec les enregistrements 

25 sismiques ; 

f) on applique une methode d'inversion cinematique avant sommation en utilisant 
les temps de trajet multi-deports obtenus pour l'ensemble des collections multi-deports 
selectionnees, afin de determiner la geometrie et la vitesse de la couche consideree ; et 
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g) on itere n fois (n>0) les etapes a) a f) en considerant a chaque iteration le 
modele de Vitesse obtenu a l'iteration prexedente comme modele de reference pour definir 
la dite vitesse de reference. 

2) Methode selon la revendication 1, caracterisee en ce que I'on realise n fois 
l'etape g) avec n>l dans le cas ou la gamme de vitesse choisie ne s'avere pas assez 
precise a Tissue de l'iteration precedente ou des etapes a) a f). 

3) Methode selon 1'une des revendications precedentes, caracterisee en ce que, dans 
le le cas ou la distribution des vitesses varie fortement lateralement et/ou dans le cas oil 
Ton ne dispose pas d'une connaissance a priori suffisamment precise sur la distribution de 
vitesse dans la couche considered, on realise l'etape g) sur des gammes de deports de plus 
en plus larges au fur et a mesure des iterations. 

4) Methode selon l'une des revendications precedentes, caracterisee en ce que, dans 
le le cas ou la distribution des vitesses varie fortement lateralement et/ou dans le cas ou 
Ton ne dispose pas d'une connaissance a priori suffisamment precise sur la distribution de 
vitesse dans la couche consideree, on realise l'etape g) sur des grilles de collections multi- 
deports de plus en plus fines au fur et a mesure des iterations. 

5) Methode selon l'une des revendications precedentes, caracterisee en ce que, a 
l'etape e), dans le cas ou la geomdtrie de l'interface recherchee est generatrice de 
triplications, on considere les branches de triplications dans les collections multi-deports 
independamment les unes des autres. 

6) Methode selon la revendication 5, caracterisee en ce que I'on realise l'etape e) en 
utilisant des outils de trace de rayons et d'inversion permettant de prendre en cornpte des 
arrivees multi-valuees. 

7) Methode selon l'une des revendications precedentes, caracterisee en ce que I'on 
construit a l'etape a) une collection iso-deport a deport nu] ou a deport proche. 

8) Methode selon l'une des revendications precedentes, caracterisee en ce que Ton 
applique a l'etape c) une technique d'inversion cinematique a vitesse fixee telle qu'une 
migration de carte. 
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9) Methode selon Tune des revendications prec6dentes, caracteris6e en ce que ]'on 
calcule a l'etape d) 1'information cinematique en tra9ant des rayons rnulti -deports sur 
1' interface de chaque couple interface-vitesse. 

10) Methode selon Tune des revendications precedentes, caracterisee en ee que Ton 
realise l'etape f) en appliquant une methode d r in version cinematique avant sommation telle 
qu'une tomographic de temps de trajet avant sommation. 
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